A model of inversion of DNA charge by a positive polymer: fractionization of the 

polymer charge 
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Charge inversion of a DNA double helix by an oppositely charged flexible polyelectrolyte (PE) is 
considered. We assume that, in the neutral state of the DNA-PE complex, each of the DNA charges 
is locally compensated by a PE charge. When an additional PE molecule is adsorbed by DNA, its 
charge gets fractionized into monomer charges of defects (tails and arches) on the background of the 
perfectly neutralized DNA. These charges spread all over the DNA eliminating the self-energy of PE. 
This fractionization mechanism leads to a substantial inversion of the DNA charge, a phenomenon 
which is widely used for gene delivery. 



Inversion of the negative charge of a DNA double hehx 
by its complexation with a positive polyelectrolyte (PE) 
is used for the gene delivery. The positive charge of DNA- 
PE complex facilitates DNA contact with a typically neg- 
ative cell membrane making penetration into the cell hun- 
dreds times more likelyEJ. Charge inversion of DNA-|PE 
complexes was confirmed recently by electrophoresisB. If 
at a given concentration of long DNA helices one in- 
creases concentration of shorter PE molecules, at some 
critical point the electrophoretic mobility of a DNA-PE 
complex changes sign from negative to positive. Intu- 
itively, one can think that when a PE completely neu- 
tralizes a DNA double helix new molecules of PE do not 
attach to DNA. Indeed, the Poisson-Boltzmann approx- 
imation for description of screening of a DNA by any 
counterions including PE does not lead to charge inver- 
sion. Counterintuitive phenomenon of charge inversion 
of a macroion by oppositely charged PE has attracted 
significant attentionB~t3. It can be explained if one takes 
into account that the surface potential of already neu- 
tralized DNA is locally affected by a new approaching 
PE molecule, or in other words, .taking into account cor- 
relations between PE moleculeaiJ'Ej. Due to repulsive 
interaction between PE molecules a new PE molecule 
pushes aside already adsorbed on DNA surface molecules 
and creates on the surface an oppositely charged im- 
age of itself. The image attract the new PE molecule 
leading to charge inversion. This phenomenon is sim- 
ilar to attraction of a charge to a neutral metal. For 
quantitative consideration charges of DNA are often as- 
sumed tei-^e smeared and to form uniformly charged 
cylinderEI~Ej. This approach ignores interference between 
chemical structure of DNA surface and of PE and clearly 
is not fully satisfactory. In this paper, we consider ef- 
fects of discreteness and configuration of — e charges of 
the DNA double helix. In this case, we suggest an expla- 
nation of charge inversion based on "fractionization" of 
charge of PE molecules. It turns out to be even simpler 
and more visual than for smeared charges of DNA. 

Negative elementary charges of DNA phosphates are 
situated along the two spirals at the exterior of both he- 
lices. When unfolded, each spiral is an one-dimensional 



lattice of such charges, with the lattice constant a=6.7A. 
Let us consider a toy model of a PE as a freely jointed 
chain of Z small +e monomers. The elastic energy cost 
for bending the PE is neglected in this model, so that one 
can concentrate on the electrostatic aspect of the prob- 
lem. To maximize the role of discreteness of DNA charge 
we assume that the PE bond length b is exactly equal to 
the distance a between negative charges of a spiral. (The 
case when these lengths are different is discussed in the 
end of the paper). We assume that minimal distance, 
d, between a PE charge and a charge of DNA is much 
smaller than a. Then PE molecules can attach to a DNA 
charge spiral in such a way that every charge of a spiral 
is locally compensated by a PE charge and, therefore, 
DNA is completely neutralized. The case of Z = 3 is 
shown in Fig. ^a. The neutralization is so perfect that it 
is difhcult to imagine how another PE molecule can be 
attached to DNA. 
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FIG. 1. The origin of charge fractionization in PE adsorp- 
tion, a) One of spirals of negative charges of DNA (empty cir- 
cles) is completely neutralized by positive PE molecules with 
Z — 3 (their monomers are shown by solid circles). A new 
PE molecule is approaching DNA. b) The new PE molecule 
is "digested" by DNA. Its charge is split in +e charges of Z 
defects. They are tails and an arch (center). 
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In this paper, we would like to discuss the fraction- 
ization mechanism which brings an additional PE to the 
neutralized DNA and leads to charge inversion. Fig. |l| 
shows how this mechanism works for the case oi Z ~ 3. 
When a new PE comes to the DNA double helix which 
is already neutralized by PE, it creates a place for it- 
self or, in other words, the oppositely charged image in 
the following way. Let us choose Z already-absorbed PE 
molecules, which are situated far from each other. In 
each of them we detach one PE monomer from DNA 
surface. This leads to formation of positive defects (tails 
and arches) and negative vacancies on DNA. To create 
an image for a new PE let us shift adsorbed PE molecules 
along DNA in such a way that all Z vacancies join to- 
gether and form a large vacancy of a length Z. A new PE 
molecule is accommodated in this vacancy. As a result of 
consumption of this molecule Z defect +e charges appear 
on the top of completely neutralized spiral (see Fig. |l|b). 

This effectively looks as cutting of the new PE molecule 
into Z individual monomers and spreading them out 
along the spiral. In other words, charge inversion of DNA 
happens by fractionization of the PE molecule charge. Of 
course, none of the chemical bonds is really cut, and this 
phenomenon is solely due to the correlated distribution of 
PE molecules, which avoid each other at the DNA spiral. 
In this sense, fractionization we are talking about is pjni- 
lar to what happens ia fractional quantum Hall effectllj or 
in the polyacetilenetll, where many-electron correlations 
result in fractionization of the electron charge. 

Fractionization is driven by elimination of the self- 
energy of free PE in solution. By the self-energy we mean 
the energy of repulsive interactions of Z positive charges 
of the PE molecule in extended conformation which it 
has in the solution. In the fractionized state, charges of 
monomers are far from each other and practically do not 
interact, so that the PE self-energy is eliminated and, 
therefore, gained. 

Let us now calculate the net inverted charge using this 
fractionization mechanism. We denote the linear charge 
density of the inverted (positive) net charge of the dou- 
ble helix DNA by 77* . The chemical potential of the PE 
absorbed at the spiral is 



Us ZkBTlii{ri* /rio) + Zei>{{)) 
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The first term in the right hand side of Eq. (|^) is the 
chemical potential of the one-dimensional gas of defects 
(— ?7o — 0.6e/A is the bare charge density of DNA). We 
used expression for the chemical potential of an ideal gas 
because the Coulomb interaction between defects at the 
a distance of a few a is much smaller than ksT (a Ib, 
where Ib = /DkBT ~ tA is the Bjerrum length.) The 
second term in the right hand side of Eq. (|l|) is the re- 
pulsion energy of the new PE from the inverted charge of 
the DNA. In this term, V'(O) is the averaged surface po- 
tential of the DNA helix. We assume in this paper that 
the net charge of DNA is screened by a monovalent salt 
at the screening length rg, which is much larger than a. 



Then V'(O) can be calculated as the surface potential of a 
cylinder with radius of DNA helix R and linear density 
of charge ry* 
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To find 77* in the equilibrium state, one has to equate the 
chemical potential of the absorbed PE molecules with 
that of the free PE in the solution. The later one can be 
calculated as following. Due to the repulsive Coulomb 
interaction between monomers, a free PE in the solution 
has an extended shape to minimize its energy. There- 
fore, the chemical potential of a free PE in solution can 
be written as the self-energy of a rigid rod with the length 
Na and the linear charge density e/a 



/io = {Ze^lDa)\n{Cla) , 



(3) 



where C = min(rs,Za) and D is the dielectric constant 
of water. We have assumed that the concentration of PE 
in the solution is large enough so that its translational 
entropy can be neglected. In this sense, we are calcu- 
lating the maximum possible charge inversion. If the PE 
molecule is long {Z ^ 1) this limit is reached at a concen- 
tration of PE which is exponentially small {'^ exp(— Z)). 

Equating the chemical potentials of Eqs. (||) and (|^), 
one has 

7/.(0) = {e/Da) \n{C/a) + (fcsT/e) ln(77o/?7*) • (4) 

One can interpret the right hand side as a "correlation" 
voltage (provided by the total free energy gain in frac- 
tionization of PE charge) that (over-) charges the DNA 
to the potential ipijS). 

To the first approximation, one can neglect the en- 
tropic term in the right hand side of Eq. (^ and easily 
get a solution for the net charge density 
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2a ln[(r« + R)/R] 
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Now one can check that this solution is consistent with 
the assumption that the entropic term can be neglected 
by substituting it back into Eq. (^. Of course, r]* is 
positive indicating that the bare DNA charge is inverted. 
Knowing 77* and using \r]o\ = 0.6e/A~ 3.9e/a the charge 
inversion ratio can be calculated 



= 0.13 



\n{C/a) 



ln[(r, + R)/R] 



(6) 



For DNA R = lOA and a = 6.7A, so that at > 
lOAthe ratio of logarithms can be only slightly larger 
than unity. Thus, the charge inversion ratio created by 
fractionization is limited by 20%. Up to such point we 
indeed can neglect Coulomb interactions between defects 
in the chemical potential of the gas of defects (the first 
term in the right hand side of Eq. (|l|) . 

We emphasize that it would be incorrect to con- 
clude that fractionization discussed above is a strictly 
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one-dimensional phenomenon, similarly to the, well 
known cases— pf one-dimensional density wavetj and 
polyacetilenell3. It is easy to see that our fractionization 
mechanism applies equally well for a two-dimensional 
surface with discrete charges which form a square lattice 
with the same lattice constant as the PE bond length c. 
Indeed, one can view Fig. |l| as a cross-section of this 
lattice and all previous arguments about the role of tails 
and arches can be carried over to this case. There are, 
however, small modifications of the analytic formulae for 
charge inversion. Defects with +e charges form now a 
two-dimensional gas with concentration cr*/e, where a* is 
the net positive surface charge density playing the roles of 
77*. The chemical potential of this gas is fcBTln(a^(T*/e). 
The surface potential is now ^jj{0) = 2na*rs/ D. The bal- 
ance of the chemical potential of PE molecules adsorbed 
at the surface with that of a free PE in the solution now 
reads 

27rcr*r,/D = {e/Da) \n{C/a) + (fcijT/e) ln(e/aV*) (7) 

and the solution, for a ~ Ib, within a numerical factor, 
is 

a* ~{e/ars)/Hrs/a) . (8) 

One can see that, for a ~ Z^, in the free energy gained 
by fractionization of the PE molecule charge, the entropy 
contribution is comparable to the self-energy, in contrary 
with the one-dimensional case, where the entropic term 
can be neglected. This is obviously due to a higher num- 
ber of degrees of freedom which a two-dimensional sur- 
face provides to the gas of defects. If Ts ^ a, the frac- 
tionization induced charge inversion ratio for the two- 
dimensional case is smaller than for DNA: 



An important role of elimination of the self-energy for 
adsorption of a flexible PE on an oppositely unifornily 
charged surface was previously emphasized in Rcfs. 

Until now we considered adsorption of linear charged 
molecules (PE) both on one- and two-dimensional lat- 
tices of the background charge. It is interesting to note 
that the fractionization mechanism works for molecules 
of other shapes, too. Let us, for example, consider den- 
drimers (star-like branching molecules with a large num- 
ber of monovalent positive charges on their peripherAti, 
which were also shown to invert the charge of DNAMI. 
Dendrimers with charges ^=4, 8 can easily compensate a 
compact group of nearest Z charges of both DNA helices. 
If a DNA double helix is totally covered and neutralized 
by such dendrimers an additional dendrimer can still be 
adsorbed on DNA because one monomer of each of Z al- 
ready adsorbed dendrimers can be raised above the DNA 
surface. As in the case of linear molecules, this leads to 
fractionization of the dendrimer charge and to the gain 
of its self-energy. 



Returning to DNA-PE complexes we would like to re- 
mind the reader about additional mechanisms of charge 
inversion beside the fractionization mechanism. So far, to 
clearly demonstrate the role of the fractionization mech- 
anism in charge inversion, we worked with the case when 
distance between charges of PE, 5, is equal to the distance 
between charges of an unfolded DNA spiral, a. One can 
show that if b varies in the vicinity of a, the point h = a 
is the local minimum of the charge inversion ratio. Away 
from b = a point, interaction of a long PE molecule with 
a spiral of DNA charges can be calculated neglecting dis- 
creteness of PE and DNA charges, i. e. assuming that 
the charge is uniformly distributed along both the PE 
molecule and the DNA spiral. 

Let us first assume that 6 < a so that the PE molecule 
has larger linear charge density than the DNA spiral. 
Then PE molecules repel each other and form on DNA 
spiral a strongly correlated liquid where PE molecules al- 
ternate with vacant places. This liquid reminds a Wigner 
crystal. A new PE molecule pushes aside previously ad- 
sorbed ones or, in other words, attracts vacancies. This 
provides another mechanism of creation of image of an 
approaching PE molecule in addition to the defect forma- 
tion mechanism. Therefore, the negative chemical poten- 
tial of PE adsorbed on the spiral becomes larger in the 
absolute value and charge inversion increases. This is 
the same mechanism of Wigne r-crystal-like correlations 
which was studied in Ref. |9[ p^ , |l^ . In the opposite case, 
6 > a, when PE has a smaller linear charge density than 
a DNA spiral, more than one layer of PE is adsorbed 
on DNA to neutralize it. Polarization of the incomplete 
second layer by a new PE molecule results again in an 
additional to defect formation mechanism of attraction 
of this molecule to a neutralized DNAEEl. This leads to 
larger charge inversion than in 6 = a case. 

Changing flexibility of a PE we can separate the role 
of fractionization. For example, for the absolutely rigid 
PE molecules defects can not appear and fractionization 
mechanism does not work. As a result at 6 = a, when 
the layer of PE neutralizing DNA is completely incom- 
pressible charge inversion vanishes (see a two-dimensional 
analog of this problem in Ref. |l2|). In a flexible PE, the 
fractionization mechanism adds charge inversion weakly 
dependent on 6, making the minimal value of the charge 
inversion ratio at 6 = a flnite. 

Wigner-crystal-like correlations play an important role 
in the discussed above case of DNA charge inversion by 
dendrimers, too. This happens when we deal with high 
generations of dendrimers which have a very large charge 
such as 32e or 64e. Because of the three-dimensional ar- 
chitecture of their chemical bonds these molecules can 
not expand enough so that each charge of them reaches 
an opposite charge of DNA and compensates it. In other 
words, when projected to a DNA double helix, these 
high generation dendrimers have much larger linear den- 
sity of charge than the double helix itself. Thus, large 
segments of the helix between adsorbed dendrimers re- 
main negatively charged, and form a Wigner-Seitz cell 
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around each dendrimer. This is how with growing charge 
of dendrimers the fractionization mechanism gets re- 
placed by the mechanism of Wigner-crystal-hke correla- 
tions. Qualitative difference between low and high gen- 
erations of depdrimers has been clearly demonstrated 
experimentallylla. 

Let us return to complexation of a DNA double he- 
lix with PE molecules with the matching bond length, 
b = a, and discuss another mechanism of charge inver- 
sion, which is related to the discreteness of DNA charge 
and further increases the positive charge of DNA-PE 
complex. Let us consider a monomer tail of PE and ex- 
plore whether some energy can be gained if the positive 
charge of this monomer moves down to the plane of DNA 
charges, approaches already neutralized negative charge 
of the DNA and shares it with the end monomer of the 
neighboring PE molecule in a way shown in Fig. ^. If 
these two end monomers may sit on opposite sides of the 
negative charge of DNA, the additional energy e^/2d can 
be gained, where d is the distance of the closest approach 
of a PE monomer and a DNA chargda. At a sufficiently 
small d this energy can be even larger than the gain per 
tail from elimination of the self-energy. In a DNA dou- 
ble helix, all the negative charges indeed are on the ridge 
above neighboring neutral atoms. Two sufficiently small 
monomers may fit into the large and small groves on both 
sides of the ridge. On the other hand if, because of stcri- 
cal limitations, they can not be in the perfect opposition 
the energy gain is smaller and can even vanish. 




FIG. 2. A view from tlie top on a spiral of negative cliarges 
of DNA (empty circles) and two PE molecules. Two positive 
end monomers share a negative charge of DNA in the perfect 
opposition. 

In this paper, we have considered three mechanisms 
of charge inversion of DNA by a PE: fractionization of 
the PE charge, Wigner-crystal-like correlations and shar- 
ing of one DNA charge by two monomers of the PE. 
We showed that depending on properties of PE they can 
work in different combinations. In conclusion we empha- 
size that all these mechanisms are due to the fact that a 
new PE molecule rearranges already adsorbed PE in such 



away that its image or correlation hole strongly attracts 
this new PE molecule. None of these mechanisms can 
be described by the Poisson-Boltzmann theory because 
this theory uses the mean-field potential which docs not 
depend on the position of the new PE molecule. Thus all 
three effects are based on correlations between different 
PE molecules. Fractionization of PE charge is the most 
visual realization of these correlations. 
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